Measuring the 13-mixing angle and the CP phase with neutrino telescopes 
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The observed excess of high-energy cosmic rays from the Galactic plane in the energy range around 
10^* eV may be naturally explained by neutron primaries generated in the photo-dissociation of 
heavy nuclei. In this scenario, lower-energy neutrons decay before reaching the Earth and produce 
a detectable flux in a 1 km"^ neutrino telescope. The initial flavor composition of the neutrino flux, 
(j>{v'e) '■ (pi^i^) • 4>ipT) = 1:0:0, offers the opportunity to perform a combined i',j/t'T appearance 
and Ve disappearance experiment. The observable flux ratio (jyipfj) / cj>{p^ -\- Ut) arriving at Earth 
depends appreciably on the 13-mixing angle i9i3 and the leptonic CP phase 5cp, thus opening a new 
experimental avenue to measure these two quantities. 



Introduction. — Neutrino physics has progressed enor- 
mously in the last decade. The discovery of neutrino os- 
cillations provides the first clear experimental signature 
for the incompleteness of the Standard Model of particle 
physics 0. While the determination of the mixing pa- 
rameters controlling the solar and atmospheric neutrino 
oscillations has already entered the precision era, it ex- 
ists currently only an upper limit for the 13-mixing angle 
i9i3 from the CHOOZ experiment, sin^ 279i3 < 0.1 
The mixing angle i?i3 characterizes how strong atmo- 
spheric and solar oscillations are coupled and therefore 
also determines the strength of CP violation effects in 
neutrino oscillations. Among the three possible phases 
in the neutrino mixing matrix, only the Dirac phase (Sep 
enters neutrino oscillations. This phase is at present com- 
pletely unconstrained. Both the mixing angle 'diz and the 
phase (5cp are observable in solar and atmospheric neu- 
trino oscillation experiments only as subleading, genuine 
three-flavor effects that are masked mainly by systematic 
uncertainties Q. While there are strong experimental 
efforts to improve the measurements of i?i3 in the near 
future by dedicated experiments 0, the detection of a 
non-zero 5cy> appears unlikely for the next generation of 
facilities [5|. Thus the answer to one of the most in- 
teresting questions in neutrino physics, namely the exis- 
tence of leptonic CP violation, probably has to await the 
construction of long-baseline experiments using second- 
generation superbeams or perhaps even a neutrino fac- 
tory. 

In the following, we propose to use high-energy neu- 
trinos produced by decaying neutrons as a new probe 
to measure 'i?i3 and (5cp with neutrino telescopes. The 
potential of neutrino telescopes to measure the atmo- 
spheric mixing angle 7?i2 has been discussed recently in 
Ref. , while tests for new physics beyond standard neu- 
trino oscillations using atmospheric neutrino data have 
been examined in Ref. 0. Previously, Refs. 0,0 dis- 
cussed possibilities to measure or to constrain i?i3 and 
i5cp for the case of decaying neutrinos, while the use of 
the neutron decay channel as neutrino source was sug- 
gested in Ref. [iflj to test quantum decoherence. Neu- 
tron primaries have been invoked to explain an excess 



of high-energy cosmic rays (CRs) from two regions in 
the Galactic plane [llj, [l^l . This signal, in a limited en- 
ergy range around 10^* eV, has been observed by several 
experiments with different techniques: The AGASA col- 
laboration found a correlation of the arrival directions of 
CRs with the Galactic plane at the 4 a level ^3 ■ This 
excess, which is roughly 4% of the diffuse flux, is con- 
centrated towards the Cygnus region, with a second hot 
spot towards the Galactic Center (GC) Q. Such a sig- 
nal has been independently confirmed by the Fly's Eye 
Collaboration and by a re-analysis of the SUGAR 
data [li |. 

Complementary evidence for a cosmic accelerator in 
the Cygnus region comes from the detection of an ex- 
tended TeV 7-ray source by the HEGRA experiment 0, 
IH | . The measured photon spectrum is difficult to explain 
in terms of electromagnetic acceleration. Also, X-ray or 
radiowave emission could not be detected by CHANDRA 
or VLA ^3 J thus favoring a hadronic accelerator. Simi- 
larly, multi-TcV 7-rays from the vicinity of the GC have 
been recently detected by HESS [20| . 

Galactic neutron sources. — The excess from the 
Cygnus and GC region is seen at £^ « 10^^ eV, i.e. at en- 
ergies where charged cosmic rays still suffer large deflec- 
tions in the Galactic magnetic field so that only a neutral 
primary can produce a directional signal. Another evi- 
dence for neutrons as primaries is that the signal appears 
just at that energy where the neutron lifetime allows neu- 
trons to propagate from a distance of several kpc. 

Neutrons can be generated as secondaries cither in col- 
lisions of high-energy protons on ambient photons and 
protons, or in the photo-dissociation of heavy nuclei. In 
the first case, the flux of Vp from neutron decays would 
be negligible compared to the neutrino flux from pion 
decays. Thus one expects a neutrino flavor composition 
of : C(i>p : (?!),- = 1 : 2 : before oscillations [13 j typical 
for most sources of high-energy neutrinos. The oscilla- 
tion phenomenology and signature for such a "standard" 
GC source were already considered in (ij. In contrast, 
photo-dissociation of heavy nuclei produces a pure v^, ini- 
tial flux. Since the energy fraction transferred to the Uf, 
is typically ^ 10"'^ and only neutrons with E < 10"'^^ eV 
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can decay on galactic distances, the neutrino flux from 
photo-dissociation is limited to sub-PeV energies. More- 
over, the threshold for photo-dissociation on UV photons 
implies a lower cut-off at E ^ TeV for the energies. 

There are several arguments in favor of the domi- 
nance of heavy nuclei in the diffuse Galactic CR flux at 
E ~ 10^^ eV. First, the end of the Galactic CR spectrum 
is expected to consist of heavy nuclei, because the Galac- 
tic magnetic field confines more easily CRs with small 
rigidity. Subtracting the spectrum expected for extra- 
galactic CRs from the measured CR spectrum, Ref. jl^l 
found evidence that the transition between Galactic and 
extragalactic CRs happens around a fcwx 10^^ eV. In this 
case, the total diffuse CR fiux between (1 — 10) x 10"^'' eV 
consists only of galactic iron nuclei and extragalactic pro- 
tons. Another method to determine the transition energ 
is to study the chemical composition of the CR flux [2S 
At present, these measurements are not fully conclusive 
but point to a dominantly heavy component in the CR 
flux at least up to ~ 10^^ eV and a possible transition to 
extragalactic protons at higher energies. Such a higher 
transition energy would also case the difficult luminosity 
requirements needed for extragalactic ultra-high energy 
cosmic ray sources . Around and above the transition 
energy, the unconflned flux from Galactic point sources 
becomes visible. If the flux from these point sources con- 
sists of protons or nuclei, has to be answered experimen- 
tally for each source separately. 

In the following we use as our basic assumption that 
photo-dissociation of heavy nuclei is the origin of the de- 
caying neutrons. We assume first that other neutrino 
sources that contaminate the pure Pe initial flux can be 
neglected, but at the end we discuss how our conclusions 
change when this assumption (that can be veriflcd exper- 
imentally) is relaxed. To be specific, we use the model of 
Anchordoqui et al. in Ref. [ll|, who calculated the neu- 
trino flux from the Cygnus region which is in the field of 
view of the km^ telescope ICECUBE . These authors 
estimated an integrated 9^ flux from neutron decays of 

2 X 10~^^cm~^s^^ at £' > 1 TeV by normahzing the 
neutron flux to the 4% anisotropic component observed 
by AGASA. This flux corresponds to « 20 events (of all 
flavors) per year in ICECUBE. 

Flavor composition after oscillations. — The fluxes (pjj 
arriving at the detector are given in terms of the proba- 
bilities Pap = P{va — > j^/s) 13^ by 



formula average-out. Then we can write 



(1) 



where we have inserted 4>a = {4>e, 0,0). Since the galactic 
distances far exceed the experimentally known oscillation 
lengths even at PeV energies, the interference terms sen- 
sitive to the mass splittings Am^ 's in the usual oscillation 
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where U is the neutrino mixing matrix and greek (latin) 
letters are used as flavor (mass) indices. 

To obtain a feeling for the dependence of the fluxes on 
iJis and (5cp, we give an expansion of Pep up to second 
order in i9i3 where we use i9i2 = | and 1^23 = f , 
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As expected, the survival probability Pee (or equivalently 
(j)^) does not depend on (5cp and the unitarity relation 
'}ZfjPep = 1 holds at each order in 1^13. Moreover, the 
and Dr fluxes depend on i5cp only via the quantity 
cos (5cp . Note that the independence of Pee from 1^23 and 
(5cp, as well as the relation Pe^ = Per {^23 1^23 + 7r/2) 
(which shows up in the op pos ite signs of the cos i5cp terms 
in Eq. lO) hold exactly [26j. Though the approximate 
relations Eq. |(2J) are useful to grasp the main features 
of the dependence of the fluxes cj)^ on 'i?i3 and 6cp, in 
the following we will use the exact expressions given in 
Eq. (|2Jl. For all numerical examples, we flx the value of 
the solar mixing angle to 'di2 = 32.5° 0. 

Flavor discrimination in ICECUBE. — Let us now re- 
call briefly the flavor-discrimination possibilities in ICE- 
CUBE m- For the energies relevant here, 10^^ eV ^ 
E < 10^'^ eV, the charged-current interactions of Ve and 
i^r are in principle only distinguishable by the different 
muon content in electromagnetic and hadronic showers. 
In practice, this is an experimental challenge and we con- 
sider Ve and Vt as indistinguishable in a neutrino tele- 
scope. By contrast, in z/^ charged-current interactions 
the long range of muons ensures that the muon track 
is always visible and allows the identification of these 
events. Finally, all flavors undergo the same, indistin- 
guishable neutral-current (NC) interactions. This inter- 
action contributes however only 20% to the total cross 
section l^^. Moreover, in this case the energy of the pri- 
mary is underestimated by a factor 3-4, further suppress- 
ing the relative importance of NC interactions because of 
the steeply falling energy spectrum. In the following, we 
neglect therefore NC interaction and consider the com- 
bined De and Dt flux (j>^ + (f>^ and the Df^ flux (j>j^ as our 
two observables. 

The flux ratio R = (j>j^/{<j>e +4'?) as only observable 
does not allow the simultaneous measurement of 1^13 and 
Sep ■ For the sake of clarity, we first explore the sensitiv- 
ity of R to the value of -dis, fixing Sep = 0. In Fig. ^ 
we show the expected ratio i? as a function of 1^13 for 
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FIG. 1: Flux ratio R = / {<f)e + <l>? ) at Earth as a function 
of i?i3 for j?23 = 35° (blue, dot-dashed curve), i923 = 45° (red, 
solid curve), ■!?23 = 55° (green, dashed curve); for initial fluxes 
0e : ^/j : (^ir = 1 : : at the source and 5cp = 0. The ratio 
R — 0.5 expected for standard astrophysical sources is shown 
for comparison. 
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FIG. 2: Flux ratio R = {(t>? + <t)? ) at Earth as a function 
of i5cp for ■&ri = 10° (red, solid curve), i?i3 = 5° (green, 
dashed curve), and ■i?i3 = 2° (blue, dot-dashed curve); for 
i?23 = 45° and initial fluxes 0e : <?!)p : 0t = 1 : : at the 
source. 



three representative values of 1923 ■ This ratio varies by 
~ 50% in the interval 0° < i^is < 10° and differs in the 
extreme by a factor of three from the standard value, 
4'^ /{4'e + 4>?) = 1/2, also shown for comparison. 

If the next generation of oscillation experiments mea- 
sures or strongly constrains iJia, a neutrino telescope may 
even aim to detect leptonic CP violation. In Fig. |21 we 
show the expected ratio R — (f>j^ /{(j)^ + 4>r) a function 
of 5cp for three values of i9i3; we have chosen the best 
fit value i923 = 45°. In this case the ratio varies maxi- 
mally by about 40% in the interval 0° < ^cp < 180° and 
differs in the extreme by a factor two from the standard 
value 1/2. If we use instead i?23 = 35° (55°), the only 
change would be an overall shift of the three curves by 
Ai?« +0.1 (-0.1). 

Event rates in ICECUBE. — The excellent angular res- 
olution of 0.7° expected for ICECUBE apphes only for 



muon induced showers, while for and Vr events the 
resolution is only about 25° According to the esti- 
mate in Ref. 0, one expects roughly 1.5 atmospheric i-'^ 
background events per year at i? > 1 TeV in a window of 
1° X 1°, i.e. « 2.3 yr~^ events in a 0.7° radius around the 
Cygnus region. This number has to be compared with the 
~ Av^ signal events assuming i?23 = 45° and 1^13 = 0. A 
2 a detection of the i?^ flux is then within 1 yr capabil- 
ity of ICECUBE. Rescaling this background number to 
a cone of 25° opening angle, one expects about 2900 z/^ 
background events and 145 background showers. Here we 
used the fact that the atmospheric neutrino background 
has a flavor ratio of (pe ■ (f>fj. ■ (f>r ^ 0.05 : 1 : in the en- 
ergy range of interest, 10" eV < £; < lO" eV H^. The 
resulting statistical fluctuation of the background shower 
number is a/TV 12. Thus integrating one year the « 16 
yr~^ rate from Cygnus one expects a 1.3 a signal hint, 
or equivalently a 4.2 cr measurement in a decade. 

Obviously, the poor angular resolution for Uf, and Vt 
events is the most serious obstacle to improve this mea- 
surement. If however a future neutrino telescope would 
be able to increase the shower resolution to, say, 10°, then 
the same estimate would lead to a 3.3 a detection already 
in one year of data taking. Theoretical predictions for the 
neutron spectrum at the source could also be used to op- 
timize the detection strategy. To fit the anisotropy data 
without introducing a cutoff, the AGASA collaboration 
required in |l3l | a source spectrum with cx or steeper, 
while the spectral index of the model of Ref. ^3 is 3.1. 
The atmospheric neutrino flux falls with a similar slope: 
its spectral index is in the range 3-3.7, being steeper 
at higher energies. Thus, if the Ve spectrum would be 
truly harder than the atmospheric neutrino background, 
the signal to background ratio could be improved by an 
increase of the threshold energy. Notice also that experi- 
mentally, the energy spectrum of the signal events could 
be more easily measured using the shower events . 

What happens to our previous estimates if we add 
some contamination from "conventional" pion decay? If 
the nuclei photo-dissociation mechanism is the correct ex- 
planation for the neutron signal, realistic models as the 
one for the Cygnus region considered in would led 
to 0(10%) flux "pollution". In this shift as low 

as 0.01-0.02 is expected in the flux ratio R, well within 
the expected experimental statistical error. An acciden- 
tal pion contamination of the same order of the expected 
signal would lead to shifts of « -1-0.1 in R: the parameter 
estimate would then be challenging, but signiflcant con- 
straints on the parameter space would be still possible, 
in particular when i^-telescopes data could be combined 
with complementary information from terrestrial exper- 
iments. Finally, we want to add a remark on the case 
when neutrons are generated mainly in pp or p7 colli- 
sions. Since the normalization of the Ve flux from neutron 
decay is based on the w 4% anisotropy in the CR data, 
the number of events in ICECUBE from neutron decay 
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does not depend on the specific generation mechanism. 
However, when neutrons are produced in pp or p^ colli- 
sions, additionally a much larger flux of neutrinos from 
pion decays with : c/)^ : « 1 : 1 : 1 is expected. 
Obviously, the background for the z?i3 and Sep searches 
discussed here would therefore drastically increase, while 
the detection of these galactic point sources by neutrino 
telescopes would become much easier. A much larger flux 
and a flavor ratio (^^/((/>f + ) « 1/2 in ICECUBE 
would be a smoking gun for the dominance of the pp or 
P7 collision mechanism. Although less exciting from the 
point of view of neutrino physics, such a measurement 
would have important consequences for the astrophysical 
source diagnostics as well as for CR composition studies 
at - IQis eV. 

Summary. — It has been argued that the excess of high- 
energy cosmic rays from the Galactic Plane in the en- 
ergy range around 10^^ cV is caused by neutron pri- 
maries generated in the photo-dissociation of heavy nu- 
clei. If this model is correct, then the initial flavor 
ratio of the neutrino flux from the Cygnus region is 
(/)e : (/)^ : 0T- « 1 : : 0. Thus Nature may provide 
in a very cheap way almost pure flavor neutrino beams, 
that similarly to proposed beta-beam factories might 
help to deepen our knowledge of the neutrino mixing pa- 
rameters. In particular, we have shown that the observ- 
able ratio 4>^ /(4>^ + 4>r) of track to shower events in a 
neutrino telescope depends appreciably on the 13-mixing 
angle and the leptonic CP phase (5cp, thus opening 
a new experimental avenue to measure these quantities. 

Obviously, a better theoretical modeling of sources as 
well as more experimental studies, not only in cosmic 
rays but also in the photon channel, are highly desir- 
able. Especially worthwhile would be a confirmation of 
the anisotropy by the Auger observatory [sJl and more 
detailed chemical composition studies by the Kascade- 
Grande experiment |32l |. 
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